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ABSTRACT 

An analytical method is presented for treating the problem of a uniformly rotating, 
self-gravitating ring without a central body in Newtonian gravity. The method is based 
on an expansion about the thin ring limit, where the cross-section of the ring tends 
to a circle. The iterative scheme developed here is applied to homogeneous rings up 
to the 20th order and to polytropes with the index n = 1 up to the third order. For 
other polytropic indices no analytic solutions are obtainable, but one can apply the 
method numerically. However, it is possible to derive a simple formula relating mass 
to the integrated pressure to leading order without specifying the equation of state. 
Our results are compared with those generated by highly accurate numerical methods 
to test their accuracy. 

Key words: gravitation - methods: analytical - hydrodynamics - equation of state 
- stars: rotation. 



1 INTRODUCTION 

The problem of the self-gravitating r ing captured the in - 
terest of such r enowned scienti s ts as iKowalewskvl l|l885l ) , 
iPoincarel (|l885al lbllcl) and lDvsonl i| 1891 1 18931 )." Each of them 
tackled the problem of an axially symmetric, homogeneous 
ring in equilibrium by expanding it about the thin ring limit. 
In particular, Dyson provided a solution to fourth order in 
the parameter a — a/6, where a provides a measure for the 
radius of the cross-section of the ring and 6 the distance of 
the cross-section's centre of mass from the axis of rotation. 
An important step toward under standing rings w it h othe r 
equations of state was taken by lOstrikeij (|l964al lbl Il965l ). 
who studied polytropic rings to first order in a and found a 
complete solution to this order for an isothermal limit. 

First numer ic al res ults for homogeneous rings were 
given by IWond (|1974| ), who was not able to clarify 
the transition to spheroidal bodies that Bardeenl 
had supposed would e x ist. | Eriguchi fc Sugimotol 



Relativity (Hachisu 1986; Ansorg, Kleinwachter & Mcincl 
l2003bl : iFisch cr. Horatschck & Ansorg 20051). 



1971 



1981 



and lEriguchi fc Hachisu] (|l98a )" developed improved meth- 
ods with which they were able to study the connection 
to the Maclaurin spheroids. Returning to the p roblem 
significantly later, lAnsorg. Kleinwachter fc Meinell (|2003d ) 
achieved near-machine accuracy, which allowed them to 
study bifurcation sequences in detail and correct erroneous 
results. It was also possible to extend the problem to non- 
homogeneous rings and even to the framework of General 
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Through the use of computer algebra, we extend 
Dyson's basic idea and determine the solution to the prob- 
lem of the homogeneous ring up to the order a 20 . We also 
present an iterative method for performing a similar expan- 
sion about the thin ring limit for arbitrary equations of state 
and general results are derived, confirmin g and generalizin g 
work that had already been published bv lOstriked (Il964bh . 
The application to polytropes is considered and ordinary 
differential equations (ODEs) are found that allow for the 
determination of the mass density. A closed-form solution 
can only be found if the value of the polytropic index is 
n — 1, and such rings are considered to the order a 3 . For 
other polytropic indices, the ODEs are solved numerically so 
that results from the approximate scheme can be compared 
to highly accurate numerical results for a variety of equa- 
tions of state. The numerical solutions considered here are 
taken from a mult i-domain spectral program, muc h like th e 
one described in lAnsorg. Kleinwachter fc Meinell (|2003al ). 
but tailored to Newtonian bodies with toroidal topologies 
(see lAnsorg fc Petrofj 120051 for more information). The so- 
lutions obtained by these numerical methods are extremely 
accurate and thus provide us with a means of testing the 
accuracy of the approximate method. 
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Figure 1. A sketch providing the meaning of the coordinates 
(r, x) in relation to the cylindrical coordinates (g, z). 



2 THE APPROXIMATION SCHEME 
2.1 The Coordinates 

To describe axially symmetric rings, we introduce the polar- 
like coordinates (r, x> y>), which are related to the cylindrical 
coordinates (g, z, tp) by 



g = b — r cos x, z = rsinx, C P = V- 



(1) 



For a given value of <p, constant values of the coordinate r 
are circles centred about (g — b, z = 0) and \ measures the 
angle along any such circle. Fig. Q] provides an illustration of 
the coordinates. The surface of the ring will be described by 
a function r s (x). In our coordinates, the Laplace operator 
applied to a function / = f(r, x) reads 



f , 1 df , 1 d 2 f 



+ 



+ ■ 



dr 2 r dr r 2 d\ 2 

/. n-1 ( df sinx df 
-(6-rcosx) cosx-j; 



(2) 



dr r d\ , 

We choose the constant 6 in such a manner that the centre 
of mass of the ring's cross-section coincide with r — 0, thus 
implying 

2tt r B (x) 

/ir 2 cosxdrdx = 0, (3) 



where fx is the mass density. 
2.2 Basic Equations 

For solving the problem of a self-gravitating fluid in equilib- 
rium, we have to fulfil Laplace's equation 



V 2 C/ ut = o 
outside the fluid and Poisson's equation 
V 2 t/ in = 4ttGm 



(4) 



(5) 

inside it, where U is the gravitational potential. Addition- 
ally, we have to satisfy Euler's equation 



dv 

H— = -/xVf/in - Vp, 
at 



(6) 



where v is the velocity of a fluid element and p the pressure. 
We consider uniform rotation about the axis g — with the 
angular velocity fl and thus have the velocity field 



leading to 



V [/ in + 



v = n x x 



-Tv-ifiV =0. 



(7) 



(8) 



We introduce the pressure function 

p 

h := 



dp' 



(9) 



which is nothing other than the specific enthalpy in the case 
of constant specific entropy. Integration gives 



U in + h- -Q 2 g 2 = Vo, 



(10) 



where Vo is the constant of integration. At the surface of the 
ring, where the pressure vanishes, we have 



Us 



1 2 2 

2 n 8 



Vo. 



(11) 



2.3 Series Expansions 

The thin ring limit is approached when the ratio of the inner 
radius Qi to the outer one g tends to 1. In this limit, the 
cross-section of the ring becomes a circle. This is the start- 
ing point of the approximation. As an interesting aside, if 
one considers a ring surrounding a central body, for exam- 
ple a point mass, the cross-section of the ring can deviate 
significantly fr om a circle even if the ra dius ratio is close to 
1, see Fig. 8 in lAnsorg fc Petrofj (120051 ). 

We describe the surface of the ring up to order q by the 
Fourier series 



(12) 



(x) =<* ( i +yzy2^ k c ° s ( k x)^ + o(<T q 



where 



(13) 



There are no sine terms because of reflectional symmetry 
with respect to the equatorial plane, a symmetry neces sarily 
present for fluids in equilibrium fsee iLichten stcin 1933]). To 
the leading order, the cross-section is indeed a circle of radius 
a. We make a similar ansatz for the mass density, pressure 
function, pressure and the potential inside the ring 



H(r,x) = Mc {^2^2fJ.ik(y)cos(kx)a+o{a q ) J , (14) 

\i=0 fc=0 / 

h(r, x) = ivG^ c a 2 I ]P ]T h ik (y) cos(fcxK + o{a q ) J , 

\i=0 fc=0 / 

(15) 



P(r, x) = irGn 2 a 2 I ^ ^Pikiy) cos(kx)a + (°" 9 ) ) , 

(16) 

Um(r,x) = -TvG^ c a 2 I ^2^2u ik (y) cos(k X )a + o(a q ) J , 

(17) 



where we have introduced the dimensionless radius 

r 



V 



(i 



(18) 



The quantity fi c is chosen to be the mass density at the point 
r — and does not represent the density's maximal value, 
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although it will not differ significantly from it in general. For 
the square of the angular velocity, we write 



Q 2 = nGn c + o{a q+2 ) J 



and for Vo 



Vo = —nG/J, c a 2 Via 1 + o(a 



(19) 



(20) 



Since the potential outside the ring contains logarithmic 
terms in r it will come as no surprise that there are In a terms 
in the coefficients of our series. Like Dyson, we introduce 



A := In — — 2. 
a 



(21) 



Because lim CT ^o o"A a = for all a, this dependence will not 
pose a problem for the iteration scheme. 



2.4 Determination of the Coefficients 

We now present a method for determining (i q k, fl q +i and 
Pqk given that the previous terms in a 1 are known. 

The idea used in Dyson's approximation scheme for ho- 
mogeneous rings makes use of the Poisson integral to deter- 
mine the gravitational potential in terms of the (s till un- 
know n) function r B (x) along the axis of rotation l|Dvsonl 
1892). This is only possible since the mass density is com- 
pletely determined for homogeneous matter once the shape 
of the ring is given. In general, however, it is necessary first 
to determine \i to the desired order before being able to 
perform the integral. By taking the divergence of the Euler 
equation ([SJl, using ([SJl and expressing h as a function of n 
using the equation of state, we obtain a second order PDE 
for fi: 



4ivGn + V 2 h-2n 2 



0. 



(22) 



Expanding in terms of a and requiring that the equations be 
satisfied for each power in a and each term in the Fourier 
expansion then results in ODEs for fJ, q k(y) once an equa- 
tion of state has been specified. These functions must be 
regular at the origin and chosen such that jiioo(O) = 1 and 
jUifc(O) — for all other i and k so as to be consistent with 
the choice /i(0, \) = /j, c . For k — 0, this condition suffices 
to determine the function uniquely. For k = 1,2, ... ,q, the 
remaining constants in the solutions of the ODEs are found 
by requiring that the pressure and thus pressure function 
vanish at the surface. Demanding this for each of the coef- 
ficients in a Fourier expansion, provides q + 1 equations for 
the remaining q constants. The additional equation can be 
used to determine /3 q o- 

Although the Poisson integral is valid everywhere, for 
technical reasons, we first calculate the potential on the axis 
of symmetry only. After this, we determine the potential 
outside the ring and in particular along its surface, where 
(fTTj) must hold. 

We label the coordinates for a point on the axis 
(R,Xr) = (r,X), irom which 



Rcos(xr) 



(23) 



follows. The axis potential is 

2tt r s (x') 

!7axis( J R) = -27VG 



2wr B ( x ') 



fj,(b — r cos x') r 



■. dr d% 



y/R 2 + r' 2 - 2Rr cos ip 
2-kG J J fi(b — rcosx) ("J^) Pi(cos ip) dr d\' 



-2-k 2 G /j, c a 2 y ] 



„ V, (21-1)!! f^) 21 - 1 Ai 



21-1 \RJ 



„ a _ 2 A-^ (2Z-1)!! fa\»-i A 

\l=l / 

(24) 

where ip :— x' ~ XR an( i Pi denote the Legendre polyno- 
mials. Please note that the expansion in terms of powers of 
1/R indicated in the last line is not trivial, since there is an 
^-dependence hidden in the terms with ip. Because of re- 
flectional symmetry, there are only terms with odd powers 
in 1/R. We expand Ai with respect to a 



Ai = 2_j auCjl + °(°" 9 ) ■ 



(25) 



Using ?7axis, we can then find the potential anywhere in the 
vacuum region. To do so, we first introduce a set of axially 
symmetric solutions to Laplace's equation that vanish at 
infinity. We define 



h(Q,z) := J 



dif 



^Jb 2 + g 2 + z 2 — 2bg cos ip ' 



(26) 



which is nothing other than a multiple of the potential of a 
circular line of mass with radius b, centred around the axis. 
In (r, x, (p)-coordinates it reads 



h(r, X ) = 



COS V I 

' Xl 46 2 -46rcosx + r 2 ) 



\J 4b 2 — 4br cos x + f 2 



(27) 



where K denotes the complete elliptic integral of the first 
kind, 



K(k) 



AO 



y/l - k 2 sin 2 6 



(28) 



Because the difference of two such solutions with different 
&'s also satisfies Laplace's equation, it is clear that 

i-i 

(29) 



where 



d d d sinx d 

db db ^ dr r &%' 



(30) 



is also a solution of Laplace's equation. Next we note that 
along the axis we have 

tt(2Z - 1)!! 



h{R) 



(31) 



{21 — l)_R 2i_1 

It then follows that the potential in the vacuum region is 
U out (r, X ) = -27vG^ c a 2 P£ a 21 ' 1 <j- l A l h + o{a q ) J , (32) 
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since this expression satisfies Laplace's equation, vanishes at 
infinity and has the correct value along the axis. For calcu- 
lating the potential at the body's surface we expand 7; (r, x) 
for r < b. For example we get 



1.00 



|)( | 86 \ | [ln(8b/r)-l]cos X r | (r 
r I 2 b V b 



(33) 



After evaluating these equations for 7; (r, x) at the surface 
r — r s (x), we use (|32|) to find the coefficients 4>ik in the 
expansion 



U s (x) = -2TvGfi c a 2 [^2^2 (pik cos(kx)cr l + o(a q 

v i=0 fc=0 



(34) 



The coefficients (j>qk still depend on q unknown f3 q k (remem- 
ber that /3 q o is already known) . By comparing the coefficients 
of cos(fcx) (k = 1,2, ... ,q), equation (jTTJ provides q equa- 
tions. Together with equation ([3} to the relevant order in 
a, we can solve for fi 9 +i and the aforementioned fi q u- The 
absolute term (cos(Ox)) gives a relation between Q, q +2 and 

Vn. 



3 GENERAL RESULTS TO FIRST ORDER 

The approximation scheme described above allows us to 
draw certain conclusions even without specifying the equa- 
tion of state, thu s generali z ing res ults that were published 
for polytropes by lOstrikerl (|l964bT l. To leading order in a, 
where nothing depends on the angle x> the ring (here a 
torus) is equivalent to an infinitely long cylinder, a prob - 
lem that was s tudied by IChandrasekhar fc Fermil (|l95St ) ; 
lOstrikerl (Il964al ). Equation (J22JI now reads 



(35) 



(36) 



(37) 

the pressure vanishes, 
r s ) =0, and we thus find 

Ml) = (38) 



d 2 Id. 

-r^ + - 3- «oo + 4^oo = 0, 
djr y ay / 



since the integrated Euler equation (I10|l tells us that 



fi 2 , 2, 



must hold, i.e. 



Oo = fii = 0. 



At the surface of the ring r 



corresponding to h(r 



and 



01 



-hi 



( d/ioc 



(39) 



y = l 



By multiplying (|35p by 7r 2 /i c o 2 fey and integrating from to 
1, one finds that the mass M to leading order can be related 
to the derivative of hoo at the point y = 1 according to 



M 



47T 2 /i. 



c a 2 b J /iooydy 




2 3-1 d/ioo 
-7r i-i c a a 



Ay 



(40) 



y = l 



A particularly interesting relation involving the square of 



0.97 -•■ 



AirbP 
GM 2 



0.94 - 



0.91 - 



0.88 



homogeneous matter 

- - polytropc with n — 1 
■ ■ ■ ■ polytropc with n — 3 



0.5 0.6 0.7 0.8 0.9 1.0 

Qi/ Qo 

Figure 2. Numerical examples demonstrating how AnbP/GM 2 
tends to 1 in the thin ring limit for various equations of state, cf. 
equation (142 1 . 



the mass can be derived by considering the integral over the 
pressure 

2tt r s (x) 

P:=2*f J P r(b-r cos X )Ar a X . (41) 
o o 

To leading order, upon taking (135[1 into account, this integral 
reads 



P = 47r 3 G/i 2 a 4 6 J pooydy 



-2^G^bJ^yUy 





r, 3^ 2 4, I d/l00 2 j 

-27T Gfi c a b J Moo-^-y ay 

o 



= 87r 3 G/i 2 a 4 & / nooy / Mooy'dj/' J dy 
o Vo / 



(42) 



47r 3 G/x 2 a 4 & / pooydy 



_ GM 2 
4%b 

Numerical examples demonstrating how AiYbP/GM 2 ap- 
proaches 1 in the thin ring limit for various equations of 
state can be found in Fig. [2] 
Equation Q tells us that 



/9uMoo(l) + J filiy ay — 
o 



(43) 



holds. The terms from the expansions (|32[) and !)25|) of the 
potential in the vacuum that we need here are 



-I 



qio = 2 1 nooydy = 



Ma 
2n 2 fj, c a 3 



(44) 



and 
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«21 = 



gMa 
27r 2 /i c a 3 



where we used (143 1) . The coefficient 0,2 from the expansion of 
the square of the angular velocity follows from the coefficient 
in front of cos(x)o" in and reads 

fi 2 =aio(l + A-2 ( 3ii) + 2a 2 i. (46) 

According to (]32|) and (I33|) . the potential at the surface is 



C M 

t/ s = -^(A + 2 + o(l)) 

7T0 



(47) 



and the constant of integration Vb then follows immediately 
from the leading order of (|lip 



\ r, = [ Uout — T^ 2 ^ 2 

GM /5A + 9 
2^6 



+ g-/3 u + o(l) 



(48) 



The term g — flu appearing in the above equation can be 
treated further by considering the rotational energy T and 
potential energy W and making use of the virial identity 



= 3P + 2T + W, 



(49) 



which then implies 



= 3P + 



GM 



nb V 5_/3U "2 



2-K 3 Gfi 2 c a i b I noohooydy. 



(50) 



By restricting ourselves to the polytropic equation of state 
(see l|88[)). we can rewrite the above integral to read 



/ 



(n + l)Kn 



l/n—l 



ivGa 2 



(51) 



(n + 1) / pooydy 



(n+l)M 2 
167TV§a 4 fr 2 ' 



where the last step follows from (|42|l . Putting this expression 
into (f50)> and using (|42j) again then yields 



g-Pi 



n-l 
8 ' 



(52) 



Taking into account 1 — Qi/ Qo = 2a, which holds to leading 
order, we can use (I46p to write 



and (|48|) can be written as 

■i-KbVo _ (-,_&_ 
5GM U V Qo 



n — 5 



5 — n 
20 



+ ln!6 



In 16 



(53) 



(54) 



for polytropes in the thin ring limit. Similar equations can 
be derived for J , T, P and via the virial identity for W (see 
Ostrikcr 1964b). These equations also hold for homogeneous 
bodies (n = 0) and numerical examples demonstrating the 
behaviour (|54[) are provided in Fig. [3] 



47!% go 

5GM 




-4 -3 
ln(l - Qi/g ) 



Figure 3. Numerical ring sequences (solid lines) for homogeneous 
matter (upper curve) and polytropes with n = 5 (lower curve) are 
plotted for rings approaching the thin ring limit. The asymptotic 
behaviour as given by equation {54) is indicated by the dashed 
lines. 



4 HOMOGENEOUS RINGS 

Homogeneous matter is defined by the simple mass distri- 
bution 



and we can write 



1 if i = k = 
otherwise 



(55) 



(56) 



The integral ([9jl for the pressure function gives h = p//J> c 
which means that hit — Pik- 

For homogeneous rings, it is possible to choose 



(57) 



flio = 
without loss of generality. 

4.1 The Zeroth Order: <x° 

To leading order, the ring is truly a torus 

r s = o[l + o(l)] (58) 

and equation <[3j is thus automatically fulfilled. The poten- 
tial along the axis is that of a solid torufl which can be 
written down explicitly 



C^axis — 



87rG/i c a 3 



3R 



1 + 



An expansion in powers of 1/R then gives 

Ax = l + o(l) =^> Q10 = 1. 



(59) 



(60) 



We can extend this axis potential to the whole exterior of 
the ring via (|32|1 , and expand it on the surface (|58|l by using 
(33} in a. We get 



t/s 



-2nGfj, c a [A + 2 + o(l) 



(61) 



1 E denotes the complete elliptic integral of the second kind, 
E(k) := J ( f \f\ - fc 2 sin 2 9d6». 
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which implies 



Plugging the expansion 



0oo = A + 2. 



(62) 



Evaluating equation confirms (|37[) and provides the re- 
lation 



Q 2 = 2« -4(A + 2). 



(63) 



This equation cannot be further evaluated until the next 
order q = 1. 

To calculate the inner structure, we have to solve equa- 
tion (|35|) . where the mass density is given by (156[) . The so- 
lution is 



hoo — -V + Ci In y + C 2 . 



(64) 



At the centre, the pressure function has to be regular, thus 
Ci = 0, and it must vanish at the surface (|38p . thus C2 = 1. 
The integrated Euler equation (|10p gives us the potential to 
this order 

Uoo = hoo + vo — , 

where both vo and £^2 are unknown, but the required com- 
bination of them is given by (|63[) and we can conclude 

Uoo = 2A + 5-y 2 , 

which indeed gives (|61|l at the surface. 

4.2 The First Order: er 1 

The surface function reads 

r s — a[l + (3ua cos x + 0(0")] 
and equation Q becomes 

2-7T 

= y rfcos^dx = a 3 [3nf3na + 0(0)]. 


This gives 

/3n = 0, (67) 



(65) 



(66) 



cf. (|43[) . which means that we have again the surface function 
(|58p . but one step further in a: 

r s = a[l + o(a)}. (68) 

Accordingly, we get the same potential as (|59|) , but with o(l) 
instead of o(a~ 1 ). The expansion in terms of powers of 1/7? 
gives 

A 1 = l + o(a) and A% = -~ + 0(0). (69) 

8 

The new coefficients are 

Q11 = and a 2 i = — — . (70) 
8 

Again we are able to calculate the potential at the body's 
surface (f68l) via (f32|) . ((33|) and the corresponding equation 
for J 2 - We get 



U s — —2nG[i c a 2 
which gives us the new coefficients 
0io = and (f>u 



A + 2+ ( - + I ) ^,,s\ 4- () [>t; 



± + 
2 8 



(71) 



(72) 



1 2 2 

2^ 8 



nGfi c a 



— + I — - fi 2 cos x ) o- + o(<t) 



(73) 

together with (I20[) and (|71[) into (|11[1 yields two new equa- 
tions when one collects the coefficients in a % cos(fcx): 



n 2 = a+ -, 
4 

^3 = 2«i. 

Using equation (|63[1 then gives 



5, 35 

2 A + y 



(74) 
(75) 

(76) 



Equation (|75|l cannot be further evaluated until the next 
order q = 2. 

For the mass and the angular momentum, we get 



and 



M = 2Ti 2 fi c a 3 [a' 1 + o{l)] (77) 



J = y/ ir 5 Gut (4A + 3) a 5 [cr~ 2 + o(<7~ 2 )]. (78) 



To leading order, fJi/go = 1 — 2(7 holds, and we can conclude 
that 



lim 

Si/go — 



AnbVo 



5GM 



In 1 



In 16, 



(79) 



as we already saw in (1541). see a lso equation (11) in 
iFischer. Horatschek fc Ansord (120051 ). 

To calculate the inner quantities, we have to find a so- 
lution to (J22J). The ansatz ([T5j) leads to the ODEs 



10 1 d/110 



2/ 



and 



d 2 /in 1 d/in 



dy 2 



+ 



2/ dy 



(80) 



(81) 



Note that the 2y term in the second equation results from 
hoo — 1 — y , which is already known. The solutions of these 
equations that are regular at the centre and obey h(r s ) = 



h 10 = 



and 



hn 



Using these results, equation (I10|) then gives 
U w = 

and 



Un = (A + l)y 



V 



(82) 

(83) 
(84) 

(85) 



After finding the inner potential and pressure, we can 
calculate the potential energy 



W 



2 



UdV 



3^-, 2 5 

-7T G/i c a 



2A 



9 



+ o(l) 



Uniformly Rotating Homogeneous and Polytropic Rings in Newtonian Gravity 7 



the rotational energy 



T = 



Q 2 dV = n 3 Gfi 2 a 5 



(x + iy-i + oia- 1 ) 



and the integral over the pressure 

P = j pdV = n 3 Gfi 2 a 5 (a' 1 + o(l)) (86) 

to first order. We see that the virial theorem ((49} is fulfilled 
at the leading order. 

For further results see Tables [BTT4B4I 



4.3 Discussion for Homogeneous Rings 

With this approximation method, we are able to calculate 
e.g. the shape, angular velocity and pressure of the ring up 
to arbitrary order in a. We have done so up to the 20th 
order. 

An important question is how good this method is. In 
Tables [T] [2] and [3] one can see how the dimensionless quan- 
tities 



M 
~M 



improve m 



p f 

accuracy 



J 
J 



1 

_ W _ 1 
" W ~ G^IqI 
with increasing 



(87) 



order for differ- 
ent radius ratios. Especially for thin rings, we get very 
accurate results. In fact, for rings with radius ratios 
Qi/ Bo ~ 0.85 we achieve a precision which is compara- 
ble with that given by the numerical m ethod described in 
lAnsorg. Kl einwachtcr & Meinel (2003a|). For larger radius 
ratios, the accuracy is thus better. As a co-product, our 
work provides an independent test of the accuracy of the 
numerical method (better than 10 -13 cf. Table [l}. 

The shape of the ring in meridional cross-section for 
various radius ratios can be found in Fig. 3] The curves to 
order q — 20 can barely be distinguished from the numerical 
ones for Qi/ g > 0.3. As one approaches the transition to 
spheroidal topologies (qi/Qo — * 0), the true curve becomes 
pointy at the inner edge and is no longer well represented by 
our Fourier series. Nevertheless, the shape of the ring is quite 
well approximated even for Qi/ g a = 0.1, as seen in Fig. [S] 
The surface function r s (x), which is a constant to leading 
order, clearly approaches the numerical one with increasing 
q. The pressure in the equatorial plane can also be seen 
to approach the numerically determined one for Qi/g a = 
0.3 in Fig. [6] It is interesting to note that the centre of 
mass does not coincide with the point of maximum pressure. 
In Fig. [7] one can get an impression of the accuracy of the 
approximation over the whole range of radius rati os and 
for va rious values of q. Despite the claims found in IWond 
|l974l ) that Dyson's perturbat ive method diverges for a > 
1/3 (see also the comments in lDvsor]|l892l ; [Bardeenlll97lf ). 
these results indicate the opposite. 



5 POLYTROPIC RINGS 

5.1 Polytropes with Arbitrary Index n 

The polytropic equation of state is 

v l+l/n 




Qi/Qo 

0.90 

0.80 

0.70 
0.50 

0.30 

0.10 
0.01 







Qo 



Figure 4. Meridional cross-sections of homogeneous rings to the 
order q = 20 for different radius ratios Qi/Qo- The q- and z-axis 
are scaled identically in such a manner that g Q has the same value 
for all the rings. The dot in each ring marks the centre of mass 
of the cross-section (g = b, z = i.e. r = 0) and the dashed 
line shows the numerical result and is indistinguishable from the 
q = 20 curve for gi/ g a ^ 0.3. 
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Figure 5. The meridional cross-section and the (dimensionless) 
surface function r s (x)/b of the homogeneous ring with radius ra- 
tio gi/g = 0.1 for different orders q compared to the numerical 
result. The surfaces are scaled such that g Q (and therefore also 
gi) has the same value to all orders. 



For large/small polytropic indices n, the equation is referred 
to as 'soft'/'stiff ' and it includes homogeneous matter as the 
limit lim n ^o K n — 1/fi — l/fi c . From now on, we shall use 
the terms 'homogeneous matter' and 'n = 0' interchange- 
ably. For polytropes, (|22[) becomes 



4nGfi + K{n + l)V 2 



- 2tr = o. 



(89) 



Instead of our coordinate y, we are now going to make use 
of a new dimensionless radial coordinate, applicable to poly- 
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Table 1. For a given radius ratio Qi/ Qo = 0.9, physical quantities to different orders in q and numerically 
determined values are compared to the values for q = 20: M 20 = 4.6299179884304816293 X 10~ 2 , tt% Q = 
3.2474683264953211610xl0~ 2 , J 20 = 7.5456215256289320669X 10" 3 , P 20 = 1.7862946528142761708X 10~ 4 , 
f 20 = 6.7988816964653749490 x 10" 4 , W 20 = -1.8956647351373578410 x 10" 3 . 
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num 
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-4 x 10~ 


14 
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14 
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14 
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14 



Table 2. For a given radius ratio Qi/q = 0.5, physical quantities to different orders in q are compared 
to the numerically determined values M num = 0.7201292, tt* um = 0.5467604, J num = 0.3247949, P nu m = 



0.04874713, T num = 0.1200820, W num = - 


0.3864053. 
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Figure 6. The pressure in the equatorial plane for a homogeneous 
ring with radius ratio Qi/ g Q = 0.3 for different orders q compared 
to the numerical result. It is interesting to note that the centre 
of mass of the cross-section does not coincide with the point of 
maximum pressure. To the order q = 18 we get b := (b— Q\)/(q — 
Qi) = 0.503 and p := (gp,max — Qi)/{Qo — Qi) = 0.480, which differ 
by less than 1% from the numerical values. 
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Figure 7. The (dimensionless) squared angular velocity f2 2 /G/i c 
as function of the radius ratio Qi/Qo for different orders q com- 
pared to the numerical result for homogeneous rings. 



and the expansion 

M = Mc 7 " (]T][> fc (f)cos(fc X K + o(a")) , (92) 

\i=0 fc=0 / 

equation (|89[) reads, cf. (1351) . 

d 2 1 d \ 4tt _ 

-7-5 + --TZ] Moo H — rMoo = 0. (93) 

dr 2 r dr ) n + 1 

This equation is sometimes referred to as one of the general- 
ized Lane-Emden equations (of the first kind) and solutions 
to it h ave been derived and studied in e.g. lGoenner fc Havasl 
(2000). No solutions other than for n = and n = 1 have 
been found for our particular parameters in closed-form and 
a discussion using symmetry tra nsformations suggests that 
they do not exist, l|Goennerl200ll ) . In the limiting case n = 0, 
this equation provides an alternative, but less transparent, 
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Table 3. For a given radius ratio Qi/Qo = 0.2, physical quantities to different orders in q are compared 
to the numerically determined values Af num = 0.9424, Q£ um = 0.9844, J num = 0.4545, P n um = 0.07865, 
f num = 0.2255, Wnum = -0.6869. 
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method to the one presented above for treating homoge- 
neous matter to lowest order. We concentrate in the next 
section on the special case n — 1. 

For other polytropic indices, the iterative method pre- 
sented here can be applied with the help of numerics. By 
describing the unknown density terms p,^ using Chebyshev 
polynomials and expanding all the quantities involved in 
terms of A, equations can be formulated for purely numeri- 
cal coefficients. The equations of the approximation scheme 
described in Section [5] must be fulfilled, whereby the ODEs 
for jlik are evaluated at collocation points of the Chebyshev 
polynomials. In general, the density functions /too are not 
analytic at f — a, meaning that high order polynomials may 
be necessary to find a good approximation of the function 
desired. We none the less chose this method, since the equa- 
tions involve integrals over the density for which one end- 
point of integration contains the unknown surface function 
r a (x), making their polynomial representation particularly 
useful. 

If one is only interested in determining a, /3n and 0,2, 
then it is not necessary to combine such numerical and alge- 
braic techniques and one can choose any numerical method 
for solving the ODEs. One begins by solving equation (|93p 
numerically for the desired polytropic index n, prescribing 
the 'initial conditions' /ioo(0) = 1 and ^/ioo| f _ = 0. For 
spherical polytropic fluids, a surface of vanishing pressure is 
known to exist only for n < 5, where the surface for n — 5 
extends out to infinity. The situation for polytropic rings is 
quite different - it seems that arbitrarily large polytropic in- 
dices are possible! Numerical solutions to (|93[) indicate that 
the density function /Ltoo indeed falls to zero for large n. The 
value of f at the first zero of the solution is a. One then 
proceeds to solve equation 



df 2 



l_d_ 

f df 



1 



An + 



47rn 

71+1 



Moo 



i - 

Mil = 



1 d/jpo 
a df 



(94) 



for /in with the condition /5ii(0) = and where jp/tii|__ 
has to be chosen so as to fulfil the centre of mass condition 
© to first order, which reads (cf. (1431) ) 



= y~ /J.iir 2 dr — n / /t l 1 fiur 2 dr 



(95) 



= a 



dr 



a/in (a) + - / floor dr. 
a 



Table 4. The values of expansion coefficients for the surface func- 
tion r s (x) and for the squared angular velocity are provided up 
to first order for different polytropic indices n. The value of a for 
n = can be found by solving J93H with n = and the conditions 
3^A*Oo|r=0 = and /ioo(0) = 1 and then locating the first zero of 
Aoo- 
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Pu 


n 2 







« 0.5642 





A + 3/4 




0.5 


0.7566 


-0.03537 


0.6371A + 0.5575 
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0.9594 


-0.07708 


0.4318(A + 1) 
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1.427 


-0.1731 


0.2169A + 0.2711 
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-3 


20 
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-2.375 


(1.362A + 7.829) x 10" 


-5 


30 


2.661 x 10 3 
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-8 
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3.337 x 10 4 


-4.875 
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50 


4.142 X 10 5 
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which now reads 



Mil (a J 



d/ioc 
df 



(96) 



and £72 is taken from (|53[) . The behaviour of these coeffi- 
cients as they depend on the polytropic index n can be found 
in Table [4] The table suggests that a — > oo and £72 — > ex- 
ponentially in n for n — > oo, which is indeed known to hold 
(|Ostri kcr 1964a, bj). The behaviour of the specific kinetic en- 
ergy of a particle in the ring, proportional to a 2 0,2 to leading 
order, will be discussed in Section [6] together with the be- 
haviour of /3n for large n. 

We provide a comparison of precise numerical values for 
various (dimensionless) physical quantities with their first 
and third order equivalents in Table [5] The dimensionless 
quantities, valid for any polytropic index n > are 



M 
M 

P 
P 



K* 



T 
T 



If 



J 
J 



G 2 p,7- 



K* 



O 
O 



Ki 



(97) 



The constant f3n can then be found using equation (|39|) . 



One can see that the accuracy of the method does not de- 
pend strongly on the polytropic index and that relative er- 
rors at third order are well below one percent for rings with 
a radius ratio of 0.9. 
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Table 5. Physical quantities to the first and third order in a are compared to the correct, numerically 
determined values for given polytropic index n and given radius ratio Q\Iq = 0.9. The polytropic indices 
n = 1.5 and n = 3 correspond to a non-relativistic and ultra-relativistic completely degenerate Fermi gas 
respectively. 
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-3174 
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num 
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5 
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-7511 
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num 
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5.2 Analytic Solution for Polytropes with n = 1 

5.2.1 The Zeroth Order: a 



We rewrite (|93[1 for n = 1, remembering that now p, = fi, 

(98) 



&r z r dr ' 



and can immediately write down the general solution 



Moo = Ci Jo(V27rf) + C 2 r (V27rf) 



(99) 



where J n is a Bessel function (of the first kind) and Y n a 
Neumann functio n (also called a Bessel function of the sec - 
ond kind), see e.g. lPrudnikov. Brvchkov fc Marichevl l)l990l l. 
The condition fi(r = 0, %) = u c tells us that Ci = 1 and 
C2 = 0. The first positive zero of Jo determines value for a 
from (|12p . We refer to the kth positive zero of the nth Bessel 
function as j„k and can then write 

a:=af/r = joi 7^2^ = 0.959... (100) 



The value for a = dy/K/G is independent of the choice 
of u c , which is not the case for other polytropic indices. 
This is due to an interesting invariance for n = 1: if U(x), 
n(x) = K^Jp(x) and v(x) are solutions to the Poisson and 
Euler equations, then so are all(x), afi(x) = K \J 'a 2 p(x) 
and *Jav(x), where a is an arbitrary scaling factor. 

5.2.2 The First Order: a 1 

The unknown quantities we have to solve for are /iio(f), 
un(f), /3io, /9n and £^2 ■ From (f2"2")l . one finds the differential 
equations 



d 2 Id, 



(101) 



and 



+ -.— ) M11 + ( 2tt- — ) fx n = --j— • (102) 



dr 2 r dr 



a dr 



Considering only solutions that vanish at the point r = 0, 
so as to maintain our choice /x(0) = /j c , we find 



uio = 



(103) 



un = C 3 Ji + -!- (V2irr Jo - Ji) , (104) 
2joi V / 



and 



where the argument of the Bessel functions is always y2nf 
unless otherwise specified. The requirement that the density 
vanish at the surface of the rings determines 

Pw = (105) 

and relates the constant C3 to the surface function 

l + 2j 2 i/3n 



C*3 = 



2joi 



(106) 



The constant /3n is determined by stipulating that the cen- 
tre of mass coincide with the point (g = 6, z = 0) as in 



o 4 - J01 

/3n "^T 



(107) 



Recalling the definition A := In — — 2, one finally obtains 



«2 = 



2JiQoi)(A + l) 
joi 



(108) 



from (fTTJ. 



5.2.3 The Second Order: a 2 

To second order, the unknown quantities that have to be 
solved for are U2o(r), p2i(f), ^22{f), foa, fhi, P22 and f23- 
The ODEs describing the mass density now read 



+ 



1 d 



dr 2 r dr 



/X20 + 2TTfl20 = 7I"^2 



2d V df f j 2d 2 dr ' 



+ 



1 d 



dr 2 r dr 



U21 + 2n 
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1 C21 _ 
r z J a dr 



and 



+ 



1 d 



dr 2 r dr 



«22 + I 2ir — — j H22 = 



1 / dun _ r d/jpo 

2d \ df f ) 2d 2 df 



(109) 



(110) 



(111) 
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The solutions vanishing at r = are 

"37Tf 2 
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/X21 = C4J1 
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(112) 
(113) 



(114) 



The constants C4 and C5 can be related to the surface 
function by requiring that ju(r s ) = hold independently for 
the coefficients in front of cos \ an( A cos 2%. The result is 



C4 = joifti 



and 



c 5 



•2 a 
-.701P22 



11 



J01 3 

64 8 4j 2 



(115) 



(116) 



Requiring the same of the coefficient in front of cos Ox gives 

1 

64" 



a _ 4 , A+l 

P20 - - — 

■7oi 



J01 



(117) 



Evaluating ((3| tells us that 

/3 2i = =► M21 = 0. (118) 
The values for the remaining constants follow from (|11[) : 

«3 = (119) 

and 



5.2-4 T/ie TViird Order: cr 3 

The third order is the final one to be presented here, but 
the iterative scheme can be applied up to arbitrary order 
assuming that one is able to solve the differential equations 
for the mass density and perform the necessary integrals. 
The ODEs that result for this order are 

d 2 1 d \ 

— + - -jz ) A*30 + 2tth 3 o = 0, 



dr 2 
d 2 



r dr 
1 d 



(121) 
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df 2 f df, 
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(122) 

(123) 



and 



2 1 d\ / 9 

+ -— )M33+ 2tt- T2 



M33 = 
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df 2 ' f df 

1 / d \ 
4a 2 V df / a V 2 df f 



r 2 d/ioo 
4a 3 df 

H22- 



(124) 



The solutions to these equations obeying the requirement 
Hqk(0) = are 



M30 = 0, 



(125) 
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(127) 
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64j 3 i 



+ 



512j 3 i 



(128) 



3(joi + 16joi(5A + 3) + 208) \ 

256^ + Cs ) J3 - 

The constants p3o, 6*6, CV and are determined by 
requiring that ^t(r = r s ) = hold independently for the 
coefficients in front of the cos Ox, cosx, cos2x and cos3x 
terms. The result is 

pso = 0, 

/A + l 3(A + 1) ^ jo 3 i 
Ge = — i ^72 J iUoi) - FTn 



2j 2 i 

• / 5 , „ 39 \ 
+ .7oi ^-A + /3 31 + — J 

14A + 9 40A + 37 



512 



32jo 



16j 3 i 



C* 7 = 
and 

C 8 



(129) 



(130) 



■hn 



/35 , 
+ I 16 A + 



•5 

■7oi 
1536 

11 



5 . _ 77 \ .3 
-A + /333 + — jjoi 



1 181 \ 119 



Equation (|3} yields 

031 



_9_ jtt 1_ 

32 + 128 jgi 



.7(U 



13, 5 
¥ A+ 16 



(131) 



1 

A V2 



-A + 9 



and (|1 ID gives 

#52 = 0, 



(132) 



(133) 
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0.25322 



Figure 8. Meridional cross-sections of polytropic rings (n = 1) 
with varying radius ratio are shown to third order (solid lines) in 
comparison to numerical results (dashed lines) for the same radius 
ratio. At the value Qi/ Q a = 0.25322 . . ., the rings reach a mass- 
shedding limit, as is evident from the numerical cross-section. 



fi 4 



o (A + ir w, >a, 
-2 12 JiUoi) + 

3oi 



+ — A + -A 

joi V § 2 



^(A+lbot 



(134) 



Joi 



11 , 17 

T A+ T 



■Mja 



and 



/ 3 33 = 7 jT-^(5A + 9) + -^ - A + 2 



64 8^ 



Joi 



(135) 



5.2.5 Physical Parameters to the Third Order 

The shape of the rings to third order that results from equa- 
tion (|12[) is compared to numerical results of the correspond- 
ing radius ratio in Fig. [8] For thin rings, the numerical and 
third order curves are indistinguishable. As the radius ratio 
is decreased, the numerical results show that the outer edge 
becomes pointier, right up to the mass-shedding limit for 
the value Qi/Q D = 0.25322 . . .. For such a ring, a fluid par- 
ticle rotating at the outer rim in the equatorial plane has a 
rotational frequency equal to the Kepler frequency, meaning 
that it is kept in balance by the gravitational and centrifugal 
forces alone - the force arising from the pressure gradient 
vanishes. The shape of the ring with the cusp that forms 
for mass-shedding configurations is not well represented by 
a small number of terms in our Fourier series. 

Using the results of the last subsection, we write down 
expressions for various physical parameters and can use 
them to verify that the virial identity (|49[) is satisfied to 



each order in a. Up to and including third order one finds 



M 



27rj5iJi(ioi 



CT 

2w 



64 



(joi + 28j 2 i + 32JQ 2 ! A - 16) Ji ( joi ) (136) 



27rjoi (A + l) Ji(joi) 



, /=4 / -W(A+l) Ji(fa)j gJJ| 

301 I cr 2 



- 2 (A+l)joiJi(joi) 



(137) 



128 (A + l) 
+ joi (96A 2 + 324A + 232) - 624A - 664] 



= V2ty . . 2 j joi 
2 a 



+ ^[-640j i (A+l) Ji (joi) 
+ 10(j ( ti-8i 2 1 + 128A + 136)]a V 



(138) 



T= ^JiOoi) 2 ioiU(A + l)j 2 i 



+ 



| - 2 (A + l) 2 Ji (joi) joi + ^ (jo 4 i (A + 1) (139) 
+ 2j 2 r(8A + 9)(2A + 3) - 112A- 132) la 



and 



w = 



•-»-;,„./,(>,)-' { - J(2A + 5) 



(140) 



- joi(2A + 5)(A + l)J 1 (joi) 

^ (jo 1 ! (2A + 5) + 4j 2 ! (16A 2 + 42A + 21) 
+ 160A + 144^ 



In the derivation of the above expressions for P and W, we 
have made use of the identity 

20 2F 3 (§, I 2, 2, |; -j O = 3j 01 2F3 (\, \; 3, 3, 7 -; - j i 



for the Gauss hypergeometric function 

pF q (01, 02, • • • ,a p ; 61,63, . . . , b q ;z) 

= (ai)fc ■ ( fl 2) fc --- (a P ) fc z fc 
' h (bi) k -(b2) k ---(b q ) k k\ 



(141) 



(142) 



with the Pochhammer bracket 



(a) h :=a(o + l)---(o + fc-l), (o) := 1. 
A proof of (|141[) can be found in Appendix [Al 
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Figure 9. The square of the dimensionless angular velocity is 
plotted versus the radius ratio for rings with polytropic index 
n = 1. 
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Figure 11. The dimensionless angular momentum divided by the 
square of the outer radius is plotted versus the radius ratio for 
rings with polytropic index n = 1. This quantity as a function of 
radius ratio tends logarithmically to infinity. 
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Figure 10. The dimensionless mass divided by the outer radius 
is plotted versus the radius ratio for rings with polytropic index 
n = 1. This quantity tends to the value M/g = 2izjoi Ji(ioi) = 
7.84 ... in the thin ring limit, which is marked by a tick. 



a case sometimes referred to as 'isothermal' because such 
an equation holds for an ideal gas at constant temperature. 
Inserting this into equation (|22p at leading order and again 
using the dimensionless coordinate f yields 

d / d 



47rr/^o 



r— In ^oo 
ar V dr 



0. 



(144) 



The solution to this equation with our normalization 
fJ-ooir = 0) = 1 reads 

Moo= (|f 2 + l)" 2 . (145) 

The density and pressure fall to zero as f — *• oo <4> r — > 
oo. Integrating over the density to calculate the normalized 
mass, one finds to leading order 



M 



floor dr = 4-7T&, 



(146) 



In order to gauge the accuracy of the expressions listed 
above, some of them are plotted to first and third order in 
comparison to numerical values in Figs[9 nTT1 The accuracy 
of the numerical values is high enough so as to render the 
corresponding curve indistinguishable from the 'correct' one 
and is plotted in its entirety, i.e. from the thin ring limit 
right up to the mass-shedding limit. The curves to first and 
third order were drawn by taking the expression for Qi, q , 
Q 2 , M and J to first and third order respectively, inserting a 
numerical value for a and then taking the appropriate com- 
bination of these numbers. One finds in all three plots that 
the third order brings a marked improvement as compared to 
the first one, but that the behaviour near the mass-shedding 
limit is not particularly well represented. 



A LIMIT OF INFINITE POLYTROPIC 
INDEX 



As n tends to infinity, the polytropic equation (|88[) shows us 
that pressure and density are proportional 

p = Kn, (143) 



which can also be read off from equation (|42[) directly, by 
making use of M = P, which is self-evident upon taking 
(|143[) into account. In Fig. 1121 the behaviour of M /4irb can 
be followed from the homogeneous case, n — 0, right up to 
the isothermal limit n — *■ oo. 

Making use of (JT45J , we find that 



lim g 

n — >oo 



0. 



(147) 



where g was defined in (|45p . It thus follows from (|52|l that 

~ (148) 



lim (f3 

n — *oo \ 



n 



as already suggested by the results of Table |U We can then 
see that the specific kinetic energy a 2 0,2 tends to infinity 
such that for fixed A 

2ira 2 Q.2 , „ 
hm = 1. (149) 

n-*oo n + 4A 

From the fact that |/3n| tends to infinity, we can con- 
clude that the range of a values for which the first order 
provides a good approximation shrinks to the point a = 0. 
In general, for a given equation of state, there exist ring so- 
lutions over an interval of radius ratios ranging from 1 in 
the thin ring limit down to a minimal value at the mass- 
shedding limit. The numerical values presented in Table [6] 
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Figure 12. The dimensionless mass divided by Arch in the thin 
ring limit is plotted versus n/(n + 1) over the whole range of 
polytropic indices n £ [0, oo). The points for n = 0, n = 1 and 
n . — > oo are known analytically and the remainder of the curve 
was generated by solving the equation for fioo numerically and 
making use of equation 1 140 I I, 



Table 6. The radius ratio of the mass-shedding ring for various 
polytropic indices n. The value Qi/ Q a = indicates the transition 
from toroidal to spheroidal topologies, which only exists for n < 
0.36. 



Bi/Qo 



0.36 




1 2 
0.25 0.37 



3 
0.44 



4 5 6 7 8 
0.49 0.53 0.56 0.58 0.60 



can conclude that ^t c oc a~ and therefore oc In a, which 
means that Q 2 and hence the velocity of a fluid element tend 
to infinity. 

Relati vistic rings, including the thin ring, limit, were 
studied in lAnsorg et all (|2003bh . lAnsorg et all (|2004r i and 
iFischer et all (|2005l ). From the perspective of General Rel- 
ativity, the Newtonian theory constitutes a good approxi- 
mation when certain conditions are fulfilled. For one thing, 
typical velocities must be small compared to the speed of 
light c and for another \U\ <C c 2 must hold. We just saw, 
however, that for rings of finite extent and mass, the ve- 
locities grow unboundedly in the thin ring limit. The same 
holds for U s as well, see l|47p . This means that the New- 
tonian theory of gravity is not appropriate to describe this 
subtle limit itself, since one cannot expect it to be a good 
approximation to General Relativity. It is remarkable that 
the approximation about the point a = is nevertheless so 
successful. 
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demonstrate that this minimal value grows for increasing n. 
The value for Qi/Q at which one reaches the mass-shedding 
limit presumably tends to 1 as n tends to infinity. This would 
imply that the isothermal thin ring limit p resented above 
and even to first order by lOstriker (|l964bh . i.e. in which 
the cross-section of the ring tends to that of a circle, is not 
unique. Analytic work including a family of isothermal thin 
ring limits and containing the mass-shedding limit will be 
presented elsewhere. 



7 CONCLUDING REMARKS 

In their work on homogeneous rings, Poincare and 
Kowalewsky, whose results disagreed to first order, both had 
made mistakes as Dyson has shown. His result to fourth or- 
der is also erroneous as we point out in the appendix. It thus 
seems particularly worthwhile to test the correctness of the 
solutions presented here. For one thing, we ensured that the 
transition condition 



viu 



VC/o, 



(150) 



is fulfilled up to the appropriate order in a. Furthermore we 
tested that the virial theorem (|49p is fulfilled for each order 
in g. 

Please note that one has to be careful in interpreting 
the results for the thin ring limit. For example, one might 
think that the squared angular velocity vanishes like a 2 In a. 
This is true for the dimensionless quantity Sl 2 /G/i c , but need 
not be true for the squared angular velocity itself. If we fix 
the 'size' b and the mass M of the ring in that limit, then 
the cross-section shrinks to a point (a — ab). With (|40[) we 
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APPENDIX A: AN IDENTITY RELATING 
HYPERGEOMETRIC TO BESSEL FUNCTIONS 

In order to prove the identity (|14ip , we prove the more gen- 
eral identity 

(§,§; 3, 3, Z;-,)- ^(§,§ ;2 , 2 ,§ : -, 

(Al) 

for an arbitrary complex number z, from which (|141|l follows 
immediately. 

We begin by using the differentiation properties of the 
hyper geometric functions, e.g. 7.2.3.47 in IPrudnikov et al] 
1 19901 ') . to write 

^ ^3 (§.f;3,3,Z;-*J z- \ 2F S (j,~;2,2,|;-z 



_d_ 

dz 



2^3 ( ^, l, 



2' 2 



(A2) 

With the integral identity 7.2.3.11, the term to be differen- 
tiated can be written as 



J (2y/UUz) zJ 1 (2 v ^t^) 



(A3) 



2 x /t 1 t a (l-*i) 2^/z(l-ti) 



dtj dt 2 



where we have made use of the identity (e.g. 7.13.1.1 in 
IPrudnikov et al.lll990l ) 



F 1 (6,- Z ) = r(6) 2 (1 - f,)/2 J 6 - 1 (2Vi). 
The above double integral yields 



(A4) 



Table Bl. Coefficients f2; for homogeneous rings up to the order 
q = 9 (i = 2,3,..., 10). 
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2^1^(1 -ti) 2^(1 -ti) 

J 

1 /■ Jo(2VrTi) . , 



dti dt 2 



(A5) 



thereby proving (|A1[I . 



APPENDIX B: FURTHER COEFFICIENTS 

In Tables [BTT4B4L coefficients for Qj, /^fe, Uik(y) and azi for 
homogeneous rings are given. 
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B2. Coefficients 0ik for homogeneous rings up to the order q 
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Table B3. Coefficients Uik(y) for homogeneous rings up to the order q = 7. 
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Table B4. Coefficients for homogeneous rings up to the order q = 8. 
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64 A ' 98304 A 
410275 \ , 25684 



1765 \3 i 1985 x 2 , 

2048 ' 1024 ' 
4356505 x , 18119345 
3538944 A + 84934656 



25 x 3 _|_ 163145 x 2 , 
384 i 1327104 ' 
62041 x , 133449991 
786432 A ~ l ~ 7644119040 



■ x 3 +-. 



■A 2 + 



A + 



47185920 

6177187 
20384317440 



10025 X 3 _ 1698785 x 2 
4096 393216 
28823195 x _ 123632507 
14155776 A 679477248 





24576 
71128549 



A - 



5308416 

12275264429 
1223059(1 Hi ill 



425 x 3 _ 5290135 x 2 _ 
12288 84934656 
82443283 x _ _gj>84148507_ 
2264924160 978447237120 



A 3 - 



6794772480 
28211731 
1141521776640 



25 X 4 , 
Wo A + 



131072 
316929175 \2 
169869312 
2116102111 
97844723712 



9525 X 3 , 

mr A + 

1001961695 
2038431744 



A + 



2195 X 4 _|_ 3052455 x 3 , 
1024 ' 524288 ' 
888829145 x 2 , 515947105 x 
169869312 ' 301989888 
13524934807 
97844723712 





20225 X 4 , 81963095 x 3 , 
98304 f 127401984 ' 
260962627 x 2 , 10908318845 x , 
382205952 ' 36691771392 ' 

1946363428441 
44030125670400 





125 X 4 , 515945 x 3 i 
24576 ~^ 42467328 ' 

1166591837 x 2 , 
101921587200 ' 

62067170651 a , 3643667826661 
12230590464000 A_l ~ 4109478395904000 



-A 2 - 



125 x 4 , 41825 x 3 i 
2359296 ~ l ~ 339738624 ' 

1276387 \2 i 692812541 x , 
12230590464 ' 17978967982080 A_r 

18307106611 
3451961852559360 



